Bifurcation at the c— 3/2 Takhtajan-Babujian point 
to the c = 1 critical lines 
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We study the S — 1 quantum spin chains with bilinear, 
biquadratic, plus bond alternation in the vicinity of the S = 
1 Takhtajan-Babujian model. Transition line between the 
Haldane and the dimer phases are determined numerically. To 
see the crossover behavior from c = 3/2 (k = 2 SU(2) WZW 
model) at the Takhtajan-Babujian point to c = 1 (k = 1 
SU(2) WZW model), we calculate the conformal anomaly c 
and scaling dimensions of the primary fields on the transition 
line. 
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I. INTRODUCTION 

Critical phenomenon is one of the most important top- 
ics in physics. For the two dimensional classical and the 
one dimensional quantum systems, the conformal invari- 
ance plays an crucial role for the understanding and the 
classification of ittl. Then the next question becomes im- 
portant: how does the crossover occur between different 
fixed points? Zamolodchikovn argued for this problem in 
general situation, and derived 'c-theorem'. This theorem 
says that if there exists a renormalization group flow from 
one (ultraviolet) fixed point to another (infrared) fixed 
point, the conformal anomaly c should decrease. With 
perturbation expansions, Zamolodchikov examiried this 
situation between minimal models (see also refjj). 

On the other hand, Affleck and Haldaneu studied quan- 
tum spin chains based on the SU(2) Wess-Zumino-Witten 
(WZW) model using non-Abclian bosonization and 0(3) 
<7 model. They argued the instability of the integrable 
Takhtajan-Babujian (TB) modeJffl which is described by 
the k = 25 SU(2) WZW modefHa. Introduction of rele- 
vant couplings by the bond alternation makes the k = 25 
system unstable and develops the excitation gap, but 
some combinations of the couplings reduce the central 
charge k (or the topological coupling of the WZW model) 
to I. In this situation, the conformal anomaly decreases 
from c — 3k/(k + 2) to c — 1, which is consistent to the 
Zamolodchikov's c-theorem. And the TB model can be 
a multi-critical pointo. 

The crossover behavior from k = 3 to k = 1 SU(2) 
WZW model was investigated by Ziman and SchulzH with 
the numerical diagonalization for the 5 = 3/2 model. 
They showed the behavior from the 5 = 3/2 TB model 
to the Hcisenberg model (with next nearest neighbor in- 
teractions). 



In this paper we discuss the other type crossover be- 
havior, that is, from k = 2 to k = 1 SU(2) WZW model. 
To see this behavior, we study the following 5 = 1 spin 
model in the vicinity of the TB point (9 = — 7r/4, 6 = 0), 
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For 6 — several authors studied numerically,0~0 and 
nowaday it is widely believed that a transition occurs 
at 9 = —it /A between the 5 = 1 Haldane (6 > — 7r/4) 
and the dimer phases (8 < —Tt/4). The ground state 
of the former phase is singlet, but for the latter phase, 
the ground state is two fold degenerate. In the Haldane 
phase, the ground state at the point 9 = arctan(l/3) 
is known rigorously.Ej This ph asa is characterized by 
the string type order parameteiilaO and the four fold 
degenerac*-|p£ the ground state for the open boundary 
conditionJlMa These are related to a hidden Z2 x Z2 
symmetry£3 In the dimer phase, the one-site transla- 
tional invariance is broken spontaneously, and neighbor- 
ing XwR spins make a singlet. Near 9 = — 7r/4, the energy 
gapD'EJ behaves as 1 9 + tt/4 | . 

With the bond alternation 5, the non-degenerate dimer 
ground state appears, and in the region — 37r/4 < 9 < 
— Xj/i—the first order phase transition occurs at S = 
0.qEJc3 There exists the transition line between the Hal- 
dane and the dimer phases. In this case, the critical prop- 
erties are of the k = 1 SU(2) WZW theory. This behavior 
has been studied for 9 = 0, 5 w 0.26 extensively-EJcil 

In this paper, we study the model ([!]) with the nu- 
merical diagonalization for finite size systems, to see the 
critical behaviors on the c = 1 critical line near the TB 
point. Our obtained phase diagram is shown in Fig. 1. 
There exist two c = 1 critical lines which separate the 
5 = 1 Haldane and two dimer phases. At the TB point, 
these two lines meet a first order phase transition line 
(9 < —tt/A, 5 = 0), on which the ground states are two- 
fold degenerate and have an excitation gap. The TB 
point is an multi-critical point and the critical behav- 
ior is described by the k = 2 SU(2) WZW model with 
c = 3/2. On the other termination point of the c = 1 
line (9,6) = (arctan(l/3), ±1), there are highly degener- 
ate ground states. 

The organization of this paper is as follows. In the next 
section, to see the instability of the c = 3/2 system, we 
consider the model based on the three component Majo- 
rana fermions as an effective one. We identify primary 
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fields of the k = 2 SU(2) WZW model with of the fermion 
model in SU(2) irreducible forms, and give some renor- 
malization group arguments. In section III, using twisted 
boundary condition method we determine the c = 1 crit- 
ical line with numerical diagonalization for finite size sys- 
tems (N — 8, 10, 12, 14, 16). To see the crossover behav- 
ior from c = 3/2 to c = 1 systems, we show the conformal 
anomaly and scaling dimensions of the primary fields on 
the critical line. Since we can only treat finite size sys- 
tems in numerical diagonalization, we analyze the finite 
size behavior of conformal anomaly in sec. II and com- 
pare it with the numerical results in subsection IIIB. The 
last section gives summary and discussion. 



II. REAL FERMIONS 

At the TB point (9 = —n/4, S = 0), the critical prop- 
erties of the system are described by the k = 2 SU(2) 
Wess-Zumino-Wittea,model.Q This has been confirmed 
by Bethe ansatzBa e3 Near this point, we have the fol- 
lowing continuum effective Hamiltonian 



H 



dr 1 

- l [.J(x).J(x): + :J(x).J(x) :] 



J dx$lf 2tl/a (x) - H J dx^ix) 
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~7T- I dx—=J(x) ■ J(x) 



(2) 



where J and J are the SU(2) currents for the left and 

right movers, the symbol means the normal ordering. 

s s z 

The operator <E> S ' s is one of the primary fields whose 
conformal weight is given by (s(s+ 1)/4, s(s+l)/4), where 
s and s are the spin of the left and the right movers, and S 
and S z are the total spin and the total third component of 
them. The last term gives a marginal correction. At the 
TB point the couplings g± and g 2 vanish, and go > (so 
that the last term is marginally irrelevant, see eq.([L2j)). 
The correspondence to the S = 1 spin model (|l) is as 
follows 



gi ocS, g 2 oc9 



(3) 



For the one-site translational invariant case in the spin 
chain (<5 = 0), g\ vanishes, and the transition between the 
two-fold degenerate dimer and the Haldane phases is in- 
duced by the third term of (||) Scaling dimension 

of the operator is x\ 1 — 1, and near the TB point 
the excitation gap behavior is AE ~ | (72 1 ~ \8 + 7r/4| as 
is the 2D Ising transition!] 

On the other hand, the k = 2 SU(2) WZW model 
can also he. described by the three component Majorana 
fermions ,Ej 
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where the fermions satisfy the anti-commutation relation 
{ip a {x),ip b {x )} = 6 a f,5(x— x'), and the SU(2) current can 
be written as 



2tt 

J a (x) = — —ie a bc^b(x)^ c (x), 

2?r - - 

J a {x) = — —te abc tp b (x)ip c (x). 



(5) 



The correlation functions of the fermions are given by 

Sab 1 



(T T [ip a (r, x)ip b (T ,x )]) 
(T T [i>a{r,x)i} b (T ,x )]) 



27T T — T + i(x — X ) 
Sgb 1 

27r r — r — i(x — x') 



(6) 



where T T is the imaginary time ordering. The operator 
product expansions (OPEs) of the current fields are 



J«(z)J b (0) = % + ^ J c (z) + • • • , 
z z z 

J a (z)J b (Q) = % + ^ J c (z) + ■■■, 



(7) 



where z = t + ix and z = t — ix. Then the primary fields 
with the conformal weight (1/2, 1/2) are written by the 
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SU(2) irreducible operators 3>f 'f as 

®i,i( x ) = (ipi(x) + itp 2 (x)) (4>i(x) + m/> 2 (£)) 

-Triip 3 (x) (4>i(x) + i$ 2 (x)) , 

—lj}\{x)^}\{x) - tf)2(x)lp2{x)) , 

-Triip 3 (x) (ipi(x) + iip2(x)) , 
${'°(x) = V2ni (Mx)Mx) ~ Mx)Mx)) , 

z £(x) = ^ (Mx)Mx) + Mx)Mx) 

+^ 3 (x)ip 3 (x)) , 

where we only show the operators with S z > (operators 
with S z < are the Hermitian conjugate of the aboiie 
ones). This fermion model was considered by TsvelikEJ 
to study the magnetic behaviox-of the 5 = 1 anisotropic 
spin chains, and Shelton et alE3 derived the model from 
the spin ladder system. 

The Hamiltonian (Q) is a simple sum of the three 
independent Majorana fermion models. By the way, 
one-component model describes the 2D Ising transition. 
Then the operator is the sum of energy densities (e). 
The order (a) and the disorder (fi) parameters are pri- 
mary fields of the c = 1/2 CFT, and one can write down 
s s z 

the primary fields 1/2 with this operators. Using the 
operator product expansionil 
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cr(z, z)fx(0, 0) 
1 



V2|z|V4 



(e OT/ V /: V' (0) + e^V^' (0)) + • • • , (9) 



(where we redefined ip a — y/2%ip a ) the SU(2) irre- 
ducible operators with the conform aJ L wejght (h, h) 
(3/16,3/16) are following two types,E 



and 



1/2,1/2 
1,0 

1/2,1/2 
0,0 

1/2,1/2 



(ciM2 + iMlO"2)M3) 
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1/2,1/2 
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1/2,1/2 
0,0 
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M1M2C3, 
Clio's, 

-^={fJ.l<J 2 + «C1M2)ct3, 
Ml^2M3- 



(10) 



(11) 



These two representations are distinguished by the parity 
of the number of order parameters. Here we choose the 
first one, so that g 2 > for 9 < — 7r/4. 

After the calculation of the OPEs among ]y 2 , 

and 3 _1 / 2 J- J, we have the following one-loop renormal- 
ization group equations, 



dlni 



rfg 2 (^) 
dbxL 



dgo(L) 
dhxL 



2- 
V3 



9i{L)g (L), 
(2-l)g 2 (L) + ^-(g 1 (L)f 



(12) 



-^y^52(i).9o(i), 
-^( 3l (i)) 2 -^( 32 (i)) 2 



(.9o(i)) 2 , 



where L is an infrared cut-off. These equations are in- 
variant under the transformation gi — > —g\ , which corre- 
sponds to 5 — > —<5 in the spin system ([!]). (On the other 
hand, these are not invariant under g 2 — > ~fQ 2 for ,gi 7^ 0. 
) In this order, Zamolodchikov's c-functiono is calculated 
as 



c( 9l (L),g 2 (L)) 



3 39 
2 ~ 32 
3\/3 



{ 9l {L)f - -{g 2 {L)f 



S 



(5iW) 2 52(i), 



(13) 



where we neglect the marginal coupling which has 
larger scaling dimension comparing two other operators. 



The RG flow is shown in Fig. 2 where we neglect go. 
It is seen from this figure, that RG flows to two different 
strong coupling regions and there is a second order phase 
transition line between them. We see that near the multi- 
critical point g\ = g 2 = Q, the second order transition line 
behaves as 



3ic~(52) 13/8 , 



(14) 



where the exponent 13/8 is the ratio (2 — x) of the two 
relevant operators. For finite L, there exists the region 
that gi{L) and g 2 (L) is small, and that the perturbation 
theory is applicable. In the leading order, the conformal 
anomaly behaves as 



c( gi (L),g 2 (L)) = ~^(g 2 (L)) 2 ---- 



(15) 



near the c = 3/2 point (a is a constant). For 51.2 7^ 0, 
the couplings gi(L) and g 2 (L) should flow to the c = 1 
fixed point or the k = 1 SU(2) WZW model. In the 
thermodynamic limit L — -> 00, the conformal anomaly 
should be discontinuous from c = 3/2 to c = 1 at the TB 
point. 



III. C = 1 MODEL APPROACH AND 
NUMERICAL RESULTS 

To investigate critical behaviors numerically, it is 
needed to determine the critical point accurately. Thus 
in this section, we firstly introduce the method to deter- 
mine the phase boundary between the Haldane and the 
dimer phases (IIIA). In subsection IIIB, we show our nu- 
merical results, and see the asymptotic behaviors of the 
scaling dimensions and conformal anomaly. 

Assuming the conformal invariance, the scaling dimen- 
sion x n is related to the energy gap of thcJipite size sys- 
tem with the periodic boundary conditionEa 



x n = -^—(E n (L) — E g (L)), 
2irv 



(16) 



where L is the system size (L = A/2), E g (L) is the 
ground state energy, and v is the sound velocity. The 
lepii|ng finite size correction of the ground state energy 



E g (L) 



TTVC 

6L 2 "' 



(17) 



where e g is the ground state energy per site for the infi- 
nite size system, c is the conformal anomaly. Using these 
relations and considering some corrections on the transi- 
tion line, we study the critical behaviors of the spin chain 
(0) numerically. 
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A. Twisted boundary method to determine the c = 1 
critical line 

In this subsection, we consider the numerical method 
to determine the critical point. Firstly let us inves- 
tigate the symmetry of the system. For the periodic 
boundary condition, the Hamiltonian ([!]) is invariant 



under two-site translation S-j 



)j+2, space inversion 
-S3 



Sj — ► Sff—j+i, an d spin reversal 5| — > —5?. Eigen- 
states have corresponding quantum numbers of the wave 
number (q = —ir — Air/N, —n — 8n/N, • • • , n), of the par- 
ity (P — ±1), and of the spin reversal (T = ±1). In a 
variational way, the ground state of theJIaldane phase 
can be described by Schwinger bosons asE3 



N-l 



,t fct 

N 



b\-b N a\) II( o 5 6 i+i- 6 i o i+i)l )' 



(18) 



while the ground state of the dimer state can be written 
as 

AT/2 

II(4i-i^-4 J --i<i) 2 lo>, (<y>o) (19) 

3=1 



JV/2-1 

J] (a^L+i-^-^+i) 2 ^), (5<0) 

where aj and are Schwinger bosons to create the 
S z = +1/2 and S z = —1/2, and |0) is the vacuum state 
for these bosons. In this boson representation, space in- 



version is P 



,t 



i! *>! & i-7+i' and s P in 



reversal is T : at <-> &t. For the periodic boundary con- 
dition, the ground states in both the Haldane and the 
dimer phases are singlet with q = 0, P = 1, and T = 1. 

To determine the transition point numerically between 
the .Haldane and the dimer phases, we use a trick in 
refs.ai23'E3 with the twisted boundary condition (TBC) 



D N+1 31 lD N+l 



-(S?±iS»), S* N+1 = Sl (20) 



This boundary condition distinguishes the above two 
states, that is, for the Haldane phase 



N-l 



(a N b\ + tf N a\) II " b i a i+i)\°)> 

3 = 1 



(21) 



which has P = T = —1, and for the dimer state (S > 0) 



N/2 

11(4 



2j-l°2j 



! |o), 



(22) 



or (S < 0) 



JV/2-1 

(aj^+fc^l) 2 J] (4/Wi- 6 2i a 2 i+ i) 2 |0), (23) 
i=i 

which have P = T = 1. These states have different 
eigenvalues for the space inversion (P) and the spin re- 
versal (T). Thus we expect that we can determine the 
phase boundary numerically by the level crossing with 
this boundary condition. 

To verify that the level crossing indicates a phase tran- 
sition, let us consider the above behavior under TBC in 
another point of view. Affleck and Haldanea argued that 
the critical properties of the considering transition point 
is of the k = 1 SU(2) WZW model. In this context, the 
phase transition is described by the following model, 



H = — I d.v 

2vr 



1 (def) 



El 

~2tt 



J dx : sin\/20 : — J dx : cos2\/20 :, 



(24) 



where n is the momentum density conjugate to 0, and to 
assure the SU(2) symmetry the couplings A-and y are 
not independent (2 A" — 2 = yo, near A = 1).E3 The dual 
field 4> is defined as d x <fi(x) = nil. Since under TBC the 
system loses the SU(2) symmetry, this model is useful 
for seeing the effect of the boundary condition. In this 
model, the operator v2 : sin \/2cf) : in the second term is 
a primary field of the k = 1 SU(2) WZW theory with the 
scaling dimension x = 1/2 and with quantum numbers 
s = 1/2, s = 1/2, S = 0. This operator is relevant and 
a second order phase transition occurs at y\ =0. The 
renormalization group equation is given by 



dyi(L) 3 3 

—— = - yii L)-- 4yi (L)y (L), 



dyo(L) 
d\nL 



1 



{ Vl {L)f - (y (L)f 



(25) 



It has been known that for a primary field : exp(zm<^ + 
in(f>) :, where integers m and n are magnetic and elec- 
tric charges in the Coulpmb gas, the effect of the twisted 
boundary conditioroocd is to shift the magnetic charge 
m as m — > m + 1/2, and n —> n. (m and n are re- 
lated to the U(l) charges a, a as a — a = \/2m, and 
a + a — y/2n.) Then the operators with the lowest 
scaling dimensions are : exp(±i0/v2) :, and the sec- 
ond term of eq. (p3) hybridizes these operators into 
0± = \/2 : cos(0/\/2± 7r/4) :. From the first order per- 
turbation calculation, size dependence of the scaling di- 
mensions of these half magnetic charge operators is given 
byEl 



x+(L) = - + 



1 , Vl (L) , y (L) 



2 

Vi(L) 



16 : 

yo(L) 

16 : 



(26) 
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respectively. yi(L) and yo(L) are the solution of eq. (25). 
Scaling dimensions are related to the excitation energy of 
the finite size system as eq. ([l6]) , thus low lying two levels 
with the twisted boundary condition ( pof ) cross linearly 
at the transition point y\ = 0. At y\ = 0, yo{L) term 
in eq. (|2^) gives a logarithmic dependence of the system 
size (yo(L) — yo/(yo\nL + 1)). Hence we verified that 
for the phase transition of the k = 1 SU(2) WZW type, 
we must observe a level crossing of two low lying energies 
under the twisted boundary conditions. 

This level -crossing is similar to the findings of Fath 
and SolyomJi-] For 6 = they studied the model (|l|) 
numerically with the general twisted boundary condition 
S^ +1 = e ±ie Sf, S Z N+1 = Sf, and they found that for 
= 7r two levels (in their notation, a and do, see Fig. 4 
in their paper) cross at the TB point. But in oup-method, 
we see different level crossing from them (in refcj, we see 
the crossing a and cq (or b)). Level crossing in ref.til 
can be related to a c = 3/2 CFT, and we can not apply 
it directly to the c = 1 case. To see this, we show the 
c = 3/2 behavior in Appendix. 



B. Numerical results 

Here we show the numerical results of N = 
8, 10, 12, 14, 16 systems with the exact diagonalization. 
We firstly determine the phase boundary between the 
Haldane a nd the dimer phases numerically (discussed in 
subsection [II A ), and next see the critical behaviors of 
the system on the c = 1 line (discussed in section II) . 

Figure 3 shows the two low-lying energies of the sub- 
space ^2 S z — with the twisted boundary condition for 
N = 12 and 8 = 0.20. These two states have different 
eigenvalues of the space inversion and the spin reversal 
operators. To compare with the periodic boundary con- 
dition case, we subtract the ground state energy E g (L). 
For 9 < 9 C , the lower state has P = T = 1, while for 
9 > 9 C the lower state has P — T = — 1. This is the 



expected behavior as is described in III A 



To see the size dependence we show the crossing points 
of N = 8, 12, 16 systems in Fig. 4. For large 8, the size 
dependence of the crossing points is small. But near the 
TB point (9 = —it/A, 8 = 0) where c = 3/2, the conver- 
gence is very slow. In Fig. 5, we show the size depen- 
dence for 8 = 0, and extrapolate data with polynomials 
of l/N. The extrapolated value is 6 = — 0.232-7T which 
deviates from the ideal value. 

To obtain the phase diagram, we extrapolate the data 
in the region 0.03 < 8 < 0.8. To calculate the scaling 
dimension and the conformal anomaly, we interpolate 
the extrapolated points with cubic spline in the region 
— 7r/4 < 9 < arctan(l/3), < 8 < 1. To compare with 
the behavior (14) of the transition line, we show the line 
8 = {9/n + 0.25) 13/8 in Fig. 4(b), but in this interpola- 
tion, we can not detect the precise exponent 13/8 of the 
critical line (|l4|). 



Near the another limit (9,8) = (arctan(l/3X-±l), there 
is the following variational calculation result .a Using the 
states ( |l8| ) and ([l9|), the transition line between the Hal- 
dane and the dimer phases is given by 
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2tan6> + 1 
'4tan0- 3' 



(27) 



Of course, this line is only valid near the point (9, 8) = 
(arctan(l/3), ±1), where states ( |l8| ) and ( |i~9| ) are exact 
degenerate ground states. In Fig. 4(a), we also show this 
line. 

We calculate the conformal anomaly on the transition 
line by the following equation 

The sound velocity v is calculated from the lowest en- 
ergy with the momentum q = Att/N and the total spin 1 
which corresponds to the SU(2) current. We determine 
the coefficients e 9 , c, and b, using N, N + 2, and N+A site 
systems. Figure 6 shows the obtained values of c(N). Re- 
flecting the c-theorem, the conformal anomaly decreases 
from c = 3/2toc=l, and increasing the system size the 
crossover region of this behavior decreases. We can see 
the asymptotic behavior of eq. (15) in Fig. 6(b). There 
exists also-Jogarithmic correction in c(L) (ground state 
energy) asoE3o 



c(L) 



1 (InL) '\ 

2+fc 8 1 ; ' 



which may be small in our numerical results (Fig. 6). 

Using eq. ([!(]), we also calculate the scaling dimensions 
corresponding to the primary fields of the k = 2 SU(2) 
WZW model in Figs. 7 and 8. These figures show the 
values of N = 16 systems with q = 0, and the total spin 
S = 0, 1, and 2 excitations. These scaling dimensions 
split due to the marginal-, corrections which depend on 
the spin of the excitationjjtj 



LAEffjL) 
2irv 



s(s + l) s(s + 1) 
k+2 + k+2 
1 [5(5+1) - s(s- 
~2 



(29) 



5(5+1)] 



InL 

In the same figure, we also show averaged scaling di- 
mensions with appropriate weights to eliminate the loga- 
rithmic correction^ Although we do not extrapolate the 
data, the averaged scaling dimensions are very close to 
the ideal value x = 3/8 and x = 1 at the TB point. We 
see that in the c = 1 region the scaling dimensions with 
s = 5 = 1/2 and 1 are £1/2,1/2 = 1/2 and x\^\ = 2 re- 
spectively. Especially the operators with s = 5 = 1 move 
to marginal ones. Comparing £1/2,1/2 an d £1,1, we see 
that Xi^i converges slowly to the expected value x = 2, 
which may come from the finite size correction of the de- 
scendent, fields of the identity operator L4I, L4Z, and 
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IV. SUMMARY AND DISCUSSION 

In this paper, we studied the S = 1 spin chains with bi- 
linear, biquadratic, plus bond-alternating interactions, to 
see the crossover behavior from the Takhtajan-Babujian 
point (9 = -tt/4, S = 0, c = 3/2) to the c = 1 
critical line. This line separates the Haldane and the 
dimer phases. Based on the k — 2 SU(2) WZW model, 
we argued the asymptotic behavior near the Takhtajan- 
Babujian point. To determine the transition point nu- 
merically, we used the twisted boundary condition on 
the basis of a c = 1 model. This method directly see the 
lower energy between two different states, which converge 
to the ground state energy of each phases in the thermo- 
dynamic limit. 

With semipJassical and topological arguments, Affleck 
and Haldanea predicted the multi-critical structure for 
spin models with bond-alternation. There should be 
25 k = 1 SU(2) WZW critical line and a first order 
phase transition line and these lines meet at the spin 
S Takhtajan-Babujian point. In the arbitrary S TB 
model, primary fields with the conformal spin h — h = 
are <f> s ,s, s — s — 0, 1/2, 1, • • • , 2S, in which operators 
with the half integer s = s have the momentum with 
q = 7T, and others have q — 0. For the one-site transla- 
tional invariant system, operators with half integer s = s 
cannot be a relevant operators in the Hamiltonian, be- 
cause these operators violate the translational invari- 
ance. Bond-alternation makes these operators into rel- 

l 2 and are 
- s s z 

>jy 2 jy 2 ShiftS tO 

1 WZW theory, and $ff 



evant one. Among these operators, $ 



1/2,1 



important for the crossover behavior. 

the primary fields in the k 
moves to the marginal operators J a J b . 

We saw the behaviors of the conformal anomaly and 
the scaling dimensions of some operators which are the 
primary fields for the k = 2 SU(2) WZW theory. This 
behavior may also be valid for bond-alternating higher 
spin cases, especially as the generalization of eq. (pT|), the 
transition line behaves as S ~ ( 52 )( 8S + 5 )/ 8S for S > 1/2 
models. 



APPENDIX A: TWISTED BOUNDARY 
CONDITIONS ON THE S = LINE NEAR THE 
TB POINT 

Fath and SolyomQ found that with the twisted bound- 
ary condition, 

Sn+i i i&N+i ~ e ±zB Si ± iSf , S%f+i = S* (Al) 

a level-crossing occurs for O = ir at the TB point on the 
line 6 = 0. Here we show this behavior, considering the 
finite size correction (up to the leading logarithmic one). 

For the twisted boundary condition, the partial 
bosonized model is useful. The critical properties of two 
independent 2D Ising models can be described by the 
Dirac fermions. This system can be identified with a 
Gaussian model. Bosonizing (||) and (|J) about the fields 
"01, tpi, "02, and 02, we have the following action 



S — So + Si, 
where the free part is given by 



(A2) 



So = 



d 2 r 1 



2^2^ 2 -' / 2. 



d 2 r 



(^3^3 +^3^3), (A3) 



where d 2 r = drdx, and we rescaled the Fermi field as 
tp — > ip = v / 27n/'. Si is the interaction part, 

Si = - 92 j d ^r)- 9 oj^l,J{r).J { r). (A4) 

A second order phase transition occurs at gi = 0. With 
the correspondence for the energy density 



-^=(iipiipi + itpzfa) = V2 : cos V2(f> 
v2 



(A5) 



and the polarization operator of the Ashkin- Teller model 
aia 2 = %/2:sin(^(/>j :, (A6) 
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we have the following operators with Sq = 



,0,0 

1/2,1/2 


= V2: 


sin 

2 V 


: o"3, 


,1.0 

1/2,1/2 


= V2: 


cos— 


: cr 3 , 




1 


^—2 : cos 


V2(t> 


*i:S 


= V2: 


sin V2(f> 




*?$ 


1 


^2^303 - 


f 2 : . 



(A7) 



In ref.0, l j 2 corresponds to b, $\'? 2 l , 2 to c , 



and 



to d . 



G 



The third component of the SU(2) current can be writ- 
ten as 



J 6 



id<p, 



J 



(A8) 



so that J 3 and J 3 are U(l) currents of the bosonic Gaus- 
sian part of eq. (|A3|), and the boundary condition (Al) 
does not affect to the fields ip3 a nd fe . For 32=90 = 0, 



the twisted boundary condition (Al) changes the opera- 
nd ■ 



tors $ x i in eq. (A7) to the following operators 
d($) = : exp 



>V2(l + £) 

Z7T 



C 2 ($) = : exp 
3 = #3^3 



n/2(-i + |^ 



(A9) 



and the corresponding excitation energy for the finite size 
system changes to 



Si(0)-B fl (O) 
E 2 (G)-E g (0) 
E 3 (e)-E g (0) 



where E g (0) is the ground state energy with the peri- 
odic boundary condition. Thus excitation levels split by 
the twisted boundary condition. The identity operator 
T (corresponding to the ground state, in refJi-1 to a) also 
changes to the operator 



2irv ( 


-) 


2tt J 


~L~ \ 


e 




2nv 

— xl > 





1(0) =: exp 
and correspondingly the energy is 

*(e)-*(o) = |(| 



(A10) 



(All) 



So far we have considered the free field theory g 2 = 
go = 0. From above equations, we see that E g (Q) and 
^2(0) are equal at = 71", and for non-zero coupling 
case, states |0 2 (0)) and |2"(0)) are hybridized by these 
couplings at $ = ir, to the operators 

1 y/2 
O e = -= : cos — </> : (= H1H2), 

= -^=:sm^<t>:{=a 1( r 2 ). (A12) 



From OPEs 



—7=-, ^O e (w,w) + 

V3 \z - w\ z 



J(z)- J(z)O e , (w,w) 



1 



1 



\/3 |z- w| 2 
1 1 



(w,w) + 
O e (w, w) 



the first order perturbation calculation of the energy is 
given by 



E e (Tr) 
_ 2tt 
~ T 

_ 2tt 
_ T 



1 1 ,n 

4 V3 V ^ 



7! ffo(L) 



(A13) 



respectively. g 2 (L) and go(^) are the solution of the fol- 
lowing renormalization group equations, 



dg2(L) 
d\nL 
dgo(L) 



92{L) 



2^3 



g 2 {L)g (L), 



3UT-T("W)"-T(»W)"- 



(A14) 



In eq. ( A13 ) ,go(i) gives the logarithmic correction at the 
TB point. Hence the energy levels a and do with $ = n, 
that is, E e (ir) and E (tt) cross linearly at the TB point 
(.92 = 0). 
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FIG. 1. Phase diagram of the S = 1 spin chains (|l|). TB 
point is 9 = — 7r/4 and 8 — (o), which is the multi-critical 
point. Lines in the region — n/4 < 9 < arctan(l/3) are 
c = 1 critical lines. • is the extrapolated numerical data 
for N = 8, 10, 12, 14, 16 systems (section III). jTihp line <5 = 
9 < — 9/4 is a first order phase transition line! 




FIG. 2. Renormalization group flow. 



N = 12, 5= 0.20 




FIG. 3. Level crossing in the system with the twisted 
boundary conditions. • is the lowest P = T = 1 energy 
E e — Eg, and o is the lowest P — T — — 1 energy E — E g for 
the N = 12 systems. We subtract the ground state energy 
with the periodic boundary condition. 
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FIG. 4. Crossing points for TV = 8 (x), TV = 12 (•), and 
TV = 16 (o) systems. +'s are the extrapolated data, and the 
solid line is the interpolated one of it. Dotted lines are (a) eq. 
© and (b) S = (0/tt + 0.25) 13/s . 
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FIG. 6. Conformal anomaly c(TV) on the c — 1 critical line 
for TV = 8 (x), TV = 10 (•), and TV = 12 (o) systems, as a 
function of (a) S c and (b) 8 C . In (b), we also show the line 
ctb — a(#+0.257r) 2 , in which numerical number ctb and a are 
determined by data at = -0.257T and -0.247T for c(TV = 12). 
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FIG. 5. Size dependence of the crossing point for 8 — as a 
function of 1/TV. The extrapolated value (•) is 8 = — 0.232-7T. 
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FIG. 7. Scaling dimensions (xf/ 2 1/2) °f ^1/2 1/2 
tained by eq. ([li]) for the TV = 16 system, x's are of 
S = and o's are of S = 1. »'s are the averaged value 
( x< i/2 1/2 + 3a;]y2 1/2)/^- Horizontal lines are 1 = 3/8 and 
x = 1/2. 
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FIG. 8. Scaling dimensions of $f'f obtained by eq. ( |li| ) 
for the N — 16 system, x 's are of S = 0, o's of S — 1, and O's 
of S = 2. »'s are the averaged value (a;?^ + 3a;} ! + 5o;i jl )/9. 
Horizontal lines are x — 1 and x = 2. 
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